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(^ . Abstract 

(N 

f— ^ ' Assuming the space dimension is not constant, but varies with the expansion of the 

^^ , universe, a Lagrangian formulation of a toy universe model is given. After a critical re- 

0^ ' view of previous works, the field equations are derived and discussed. It is shown that 

Q I this generalization of the FRW cosmology is not unique. There is a free parameter in the 

O^' theory, C, with which we can fix the dimension of space say at the Planck time. Differ- 

^ I ent possibilities for this dimension are discussed. The standard FRW model corresponds 

bX)' to the limiting case C —t +oo. Depending on the free parameter of the theory, C, the 

expansion of the model can behave differently to the standard cosmological models with 
constant dimension. This is explicitly studied in the framework of quantum cosmology. 
. , The Wheeler-De Witt equation is written down. It turns out that in our model universe, 

j^ ■ the potential of the Wheeler-DeWitt equation has different characteristics relative to the 

potential of the de Sitter minisuperspace. Using the appropriate boundary conditions 
and the semiclassical approximation, we calculate the wave function of our model uni- 
verse. In the limit of C — > +00, corresponding to the case of constant space dimension, 
our wave function has not a unique behavior. It can either leads to the Hartle-Hawking 
wave function or to a modified Linde wave function, or to a more general one, but not to 
that of Vilenkin. We also calculate the probability density in our model universe. It is 
always more than the probability density of the dc Sitter minisuperspace in 3~space as 
suggested by Vilenkin, Linde, and others. In the limit of constant space dimension, the 
probability density of our model universe approaches to Vilenkin and Linde probability 
density being exp(— 2|S'£;|), where Se is the Euclidean action. Our model universe indi- 
cates therefore that the Vilenkin wave function is not stable with respect to the variation 
of space dimension. 
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1 INTRODUCTION 

Kaluza-Klein and string theories are well known for allowing the space dimension to be other 
than 3, but still an integer and a constant ||l|, ^. This, being considered for the high energy 
limit in the universe or for the dimension of space at the Planck time, has encouraged people 
to suggest that the dimension of space in the lower energy limit, or for the actual structured 
universe, is other than 3. Recently models have been proposed in which the universe has 
(1 + 3 + re)-dimensional space with Planck scale near the weak scale, with n > 2 new sub- 
millimeter sized dimensions [y]. 

In the last years there has been attempts to identify a fractal dimension for the matter dis- 
tribution in space using either the cosmic microwave background radiation or the galaxy 
distributions Q, ^. Aside from the actual dimension of space or the matter distribution in 
it, it is interesting to study the cosmological consequences of a fractal and variable space 
dimension. All critiques of space dimensionality other than 3 rely upon cosmologically small 
scale observations [^. Therefore, one could ask about the consequences of a dynamical space 
dimension in cosmological time and space scales. A proposed way of handling such a concept 
is using the idea of decrumpling coming from the polymer physics [7-9]. The evolusion of the 
fractal dimension of a self-similar universe in the context of Newtonian gravitation is studied 
in [|10|. 

Our toy model universe consists of small cells having arbitrary space dimensions. The ef- 
fective spatial dimension of the universe in large will depend on the configuration of the 
cells. As an example, take a limited number of small three dimensional beads. Depending 
on how these beads are embedded in space they can configure to a one-dimensional string, 
2-dimensional sheet, or a three dimensional sphere. This is the picture we are familiar from 
the concept of crumpling in the polymer physics where a crumpled polymer has a dimension 
more than 1. Or take the picture of a clay which can be like a 3 dimensional sphere, or a 2 
dimensional sheet, or even a one dimensional string; a picture we are familiar from the theory 
of the fluid membranes. In this picture the universe can have any space dimension. As it 
expands the space dimension decreases continously. 

In this paper we modify and study throughly the model proposed in [11-13]. The original 
decrumpling model of the universe seems to be singularity free having two turning points for 
the space dimension |l^, |l^ . Lima et al. [|T^ criticize the way of generalizing the standard 
cosmological model to arbitrary space dimension used in [^ and propose another way of 
writing the field equations. Their model shows no upper bound for the dimension of space. 
As the universe expands the spatial dimension decreases to D = 1 [13]. 



Later on this scenario was extended to the class of multidimensional cosmological models, 
where extra factor spaces play the role of the matter fields. In this multidimensional cosmo- 
logical model an inflationary solution was found together with the prediction that the universe 
starts from a nonsingular space time [0]. Melnikov et al. has written the Wheeler-DeWitt 



equation for multidimensional cosmological model in any constant spatial dimension [15|, see 



also [16 1 and references therein. 

We review briefly the previous works on the decrumpling models and show some difficulties 
in their formulations and conclusions. To remove these difficulties, we propose a new way to 
generalize the gravitational action in constant dimension to the case of dynamical dimension. 
We use the Hawking-Ellis action for a prefect fluid in 3~space [^ and generalize it to the 
case of dynamical dimension. As we will show in details the generalization of the gravita- 
tional action to the dynamical dimension is not unique. Moreover, in contrast to the earlier 



works in [11-13], we take into account the dependence of the measure of the action on space 
dimension. The generahzation of the action, the Lagrangian, and the equations of motion 
to dynamical space dimension is then done in two ways. Studing the time evolution of the 
spatial dimension, we obtain some numerical results for the turning points in our model. 
In the second part of the paper we study the quantum cosmology of our model universe. The 
Wheeler-DeWitt equation is written and the zero points of its potential are discussed. It is 
worth noticing that the potential of our model has completely different behavior from the 
potential of the de Sitter minisuperspace in 3-space. Imposing the appropriate boundary 
condition in the limit a — > +00, and using the semiclassical approximation, we obtain the 
wave function of our model. It is then seen that in the limit of constant space dimension, 
our wave function is not well-defined. We show that it can approach to the Hartle-Hawking 
wave function or to the modified Linde wave function, but not to that of Vilenkin. We also 
estimate the probability density in our model universe. In the limit of constant spatial di- 



mension, it approaches to Vilenkin, Linde and others' proposal [18]; i.e. to the probability 
density V oc exp(25£;), or more generally exp(— 2]S'£;]), where Se is the Euclidean action of 
the classical instanton solution. 

The paper is organized as follows. In Sec. g we discuss the dimensional constraint, proposed in 
]11]. In Sec. y, we briefly review the earlier works on the decrumpling model and discuss their 



shortcomings. In Sec. H, a new way to generalize the gravitational- and the matter-action to 
the case of dynamical dimension is presented. The equation of motion and the time evolu- 
tion equation of spatial dimension is then obtained. Sec. || is devoted to the Wheeler-DeWitt 
equation of our model and the calculation of its wave function in semiclassical approximation. 
The probability density in our model universe is then estimated semiclassically. Section ^ 
concludes the discussion of our toy model universe. 

2 MOTIVATION AND CONSTRAINT OF THE MODEL 

There are observational evidences that the matter distribution in the universe, up to the 
present observed limits of 100/i~^Mpc, is a fractal or multifractal [0, |5| having a dimension 
of about 2. This should be considered as an effective dimension of space without questioning 
the three dimensionality at cosmologically small scales 0]. 



To interpret the effective space dimension D we follow the picture proposed in ]ll]-]|l^ 
Imagine the fundamental building blocks of the universe as cells being arbitrary dimensional 
and having, in each dimension, a characteristic size 6 which maybe of the order of the Planck 
length 0(10^^'^ cm) or even smaller. These "space-cells" are embedded in a V space, where V 
may be up to infinity. Therefore, the space dimension of the universe depends on how these 
fundamental cells are configured in this embedding space. The universe may have begun from 
a very crumpled state having a very high dimension D and a size 5, then have lost dimension 
through a uniform decrumpling which we see it as a uniform expansion. The expansion 
of space, being now understood as a decrumpling of cosmic space, reduce the space-time 
dimension continuously from D + 1 to the present value Dq + 1. We do not fix Dq to allow 
being other than 3. 

Here, we are interested in the correlation between the radius of our model universe and its 
dimension. In the next section, we will implement this idea in the Lagrangian formulation. 
The first major formal difficulty in formulating a space-time theory with variable space 
dimension is that the measure of the integral action is variable and therefore some part of 



integrand. The variational calculus for such a case has not been formulated yet. We are, 
however, interested in a cosmological model for the actual universe. Therefore, it is reasonable 
to accept the cosmological principle: the homogeneity and isotropy of space. This simplifies 
the matter substantially, so that it becomes possible to formulate a Lagrangian for the theory 
and write the corresponding field equations. 

2.1 Relation between the effective space dimension D{t) and characteristic 
size of the universe a{t) 

Assume the universe consists of a fixed number A^ of universal cells having a characteristic 
length 5 in each of their dimensions. The volume of the universe at the time t depends on 
the configuration of the cells. It is easily seen that |l^, |l2[ 

volD(cen) = volDo{cell)d^~^° . (1) 

Interpreting the radius of the universe, a, as the radius of gyration of crumpled "universal 
surface" Q, the volume of space can be written [0, 12]: 



a^ = iVvolD(cen) 

(2) 

or 

i-gf = e^, (3) 

where C is a universal positive constant. Its value has a strong influence on the dynamics 
of space-time, for example on the dimension of space say at the Planck time. Hence, it has 
physical and cosmological consequences and may be determined by observation. We coin 
the above relation as "dimensional constraint" which relates the "scale factor" of our model 
universe to the space dimension. Aside from any dynamics, which we will discuss in the next 
sections, it is worth discussing the dimensional constraint thoroughly. 

2.2 Discussion on the dimensional constraint 

Eq.(y) can be written as 

The time derivative of the above relation leads to 

In the above relation qq is the present scale factor of the universe corresponding to the 
space dimension Dq. The values of C and 5 can be calculated now in terms of other known 
quantities: 

[D - Dq) a 



and 



log(- 



logTpr 


1 -U D 



(7) 



where we have replaced a^ by the present value of the Hubble radius Hq ^ 



3000/i'iMpc = 
9.2503 X 10^^/i~^cm. As the values of C and 5 are not very sensitive to 0.5 < /i < 1 we take 
for simplicity h = 1. 

Note that for D -^ +oo the radius of the universe, a{t), tends to 6 (see Eq.|7|). In principle, 
6 may be greater or less than the Planck length Ipi. Assuming 5 = Ipi = 1.6160 x 10~^^cm 
and taking Dq = 3 we obtain from (y) C = 419.7. Values of 5, or the minimum radius of 
the universe, less than the Planck length correspond to C > 419.7. Assuming C < 419.7, 
we obtain for the minimum radius of the universe, 5, values bigger than the Planck length. 
Not all dimensions of space at the present or the Planck time are of interest. We are just 
interested in the experimental value Dq = 3— (5.3ib2.5) x 10^^ [19| based on the g — 2 factor of 
electron, or Dq ~ 2 as fractal dimension for matter distribution in the universe coming from 
cosmological considerations Q. At the Planck time we assume the space dimension, Dpi, to 
be about 3 or one of the values 4, 10, or 25 coming from supergravity ||2| and superstring 
theories pO| , ^ . It is reasonable to assume that the minimum radius of the universe is less 
than the Planck length. Therefore, in the following we will assume C > 419.7. Note that, 
in an expanding universe, the positivity of C gives the decreasing of the spatial dimension 
in the course of time (cf. Eq.^). Table I shows values of C and S for different values of Dpi 
assuming Dq ~ 3. Similarly, table II shows the same values for Dq = 2. 
Figs. 1-4 demonstrate changes of C or 6 vs Dpi and Dq. Figures 1-2 demonstrate the 
changes of C vs Dpi (or log Dpi) for different values of Dq in different ranges. In order to 
illustrate the behavior of 6 vs Dpi, Figs. 3 and 4 are drawn. In general, for a fixed value of 
Dq, when C increases 6 and Dpi decrease. 

TABLE I. Values of C and S for some interesting values of Dpi assuming Dq ~ 3. The reason for 
taking Dq = 3 — 5.3 x 10~^ is based on the experimental data by measuring g-factor of electron (for 
details see Ref [|l9|] ) . Dpi could be any value, here we take an arbitrary value for it, of course with 
condition Dpi > Dq. 



Dq 


Dpi 


C 


(5 (cm) 


3-5.3X 10-^ 


3 + 10-^ 


1.2588 x 10^ 


4.4747 X io-i«^^i" 


3 


4 


1.6788x10^ 


8.6158 X 10-^1^ 


3 


10 


5.9957 X 10^ 


1.4771 X 10-^*^ 


3 


25 


4.7693 X 10^ 


8.3811 X lO-'^^ 


3 


+00 


4.1970 X 10^ 


1.6160 X 10"^^ 



TABLE II. Values of C and S for different values of Dpi assuming Dq = 2. 



Dpi 


C 


(5 (cm) 
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5.5960 X 10" 


2.8231 X 10"-^* 


10 


3.4975 X 10" 


1.0447 X lO-'i^ 


25 


3.0413 X 10" 


8.4171 X 10-^'^ 


+CX3 


2.7980 X 10" 


1.6160 X 10~^^ 



3 REVIEW OF THE DECRUMPLING UNIVERSE MODEL 



The assumption of variable space dimension brings in serious difficulties in formulating field equa- 
tions. To formulate a gravitational theory based on geometry we do need a tensor calculus which 
is not defined for dynamical space dimension. It may be possible to write an action but we do not 
know how to formulate a calculus of variation in the case that the measure of the integral is itself a 



dynamical variable. Therefore, we have to look for methods to define or obtain the field equations. 
The homogeneity and isotropy of space-time, being the main characteristics of the Friedmann models, 
allow us to formulate an action and a variational method to obtain the field equations. Although this 
do not leads to a unique model for a universe having a continuously varying dimension, it leads to 
interesting results and maybe the simplest way to implement the idea of having a space with variable 
dimension. 

3.1 Action and Lagrangian 

As was mentioned before, we have no way of defining an action for a general metric implementing the 
dynamical space dimension. Therefore, let us define it first for the special FRW metric in an arbitrary 
fixed space dimension D and then try to generalize it to variable dimension. Now, take the metric in 
constant D + 1 dimension in the following form (we use natural units or high energy physics units in 
which the fundamental constants are h = c ^ ks ~ 1, G ^ Ipi = \/mpi^). 

ds^ = -N^{t)dt^ + a^{t)dllk^, (8) 

where N{t) denotes the lapse function and dY^k is the line element for a D~manifold of constant cur- 
vature k — +1, 0, —1, corresponding to the closed, fiat, and hyperbolic spacelike sections, respectively. 
The Ricci scalar is given by 

D ,2a , ,,,d,^ N^k, 2dN ^ , , 

N'^ a a a^ aN 

For simplicity, let us now take fc = in our treatments as done in Refs. |0, 13|. Substituting from 
Eq.(0) in the Einstein-Hilbert action for pure gravity, 

SG = ^Jd^'+'''>xV^R, (10) 

where k = SttG, we are led to 



a 



— d^^^ 'x {total time derivative 
2k J ^ 



N 
The Lagrangian then becomes 



DiD-l)^d^,^^^_ (U) 



L°^ = -^fe^(-f «". (12) 



2kN 'a 
Following Ref . mU , we introduce the matter Lagrangian for a prefect fluid as 



<> = -^ + ?i^, ,13) 



where 



P ■■- ^a^, (14) 

p := pNa^-^S'^, (15) 

and p and p being the energy density and pressure, respectively. Therefore, the complete Lagrangian 
is 

L^"^--^^C-ra-H-^ + '^). (16) 

Zki\ a 2 2 



This Lagrangian suffers from the fact that its dimension is not (length) ^^ as we know it from the 
D = 3 case, and depends on the space dimension. To remedy this shortcoming which is an obstacle 
for generalizing to variable dimension, we first assume that the dimension of k is constant and equal 
to (length)^""^ for Dq = 3. This is equivalent to the assumption that for k we take the usual three 
dimensional gravitational constant. We then multiply the above Lagrangian by gq^"'^, where qq 
is the present scale factor of the universe. For simplicity, we have assumed that a and ao have the 
dimension of length. Omitting the constant factor oo^" we arrive at 

/: = -^^A^(-)- + (-^ + ^), (17) 

Zk1\ a ao 2 2 

where 

P=^(-)^, (18) 

N ao 

p^pa-^N{—)^. (19) 

ao 

This Lagrangian has the dimension (length)"'* which is due to omitting of the factor ao^° for Do — 3. 
It is now in a form ready to generalize to any variable dimension. We consider the gravitational 
coupling constant k as a constant throughout this paper. The dependence of k on the constant spatial 
dimension is studied, in Ref. |2j. 

3.2 Field equations 

For D = Do, the Friedmann equations derived by the Einstein field equations are 

N-' a Do[Do - 1) 

and 



■^^ = 77775 T7 [(2 - i5o)p - Dop] , (21) 

TV^a Do[Do - 1) 

with k = 0. Note that the field equations are written in the gauge N = 0. Eqs.(|2y, |l]) are the 00- 
and non-vanishing ij-component of the Einstein field equations, respectively. There is an alternative 
Lagrangian method of deriving these field equations which can be generalized to the case of dynamical 



space dimension. Consider first the usual case of Z? = Do = const. It is easy to show that (20) is the 
equation of motion corresponding to variation of the Lagrangian (RT^ with respect to N. Moreover, 
variation of ( [l7| ) with respect to a yields 

-iDo-l)[^ + {^^){^f]^npN^ (22) 

a 2 a 

which is a combination of (|2y and pi|) . Note that in varying the Lagrangian with respect to A'^ and a the 
quantities p and p are considered as constants |ll| - |13| . Therefore, for constant space dimensions this 
Lagrangian method leads to the familiar Friedman equations derived by the Einstein field equations. 
The continuity equation for prefect fiuid is obtained from the Friedman equations (^ and |2^) or ( po| 
andU): 

|(p«^")+p|(a^'')=0. (23) 

This, in addition to an equation of state, completes the field equations for a flat FRW universe for 
any fixed space dimension Do. 

To generalize the cosmological model to spaces with variable dimension, we are in principle free to 
choose any of the above approaches. But, it is worth noticing that the tensor calculus is meaningless 
in the case of dynamical dimension, since the index of a tensor is always constant and integer number. 
Therefore, there is no way of generalizing the Einstein equations directly. In contrast, it is possible 



to write down generalized Lagrangian, such as (|17[), due to the fact that space is homogeneous and 
the dynamical part of the measure is well defined and separated. Hence, we are left with the second 
approach as the only possible way of formulating a cosmological model incorporating the dynamical 
nature of space dimension. However, we are faced now with different kind of non-uniqueness in 
formulating the Lagrangian or the field equations, which are discussed in the following sections. 

3.2.1 Original decrumpling model 



Following [y_l|, we vary the Lagrangian (|17|) with respect to a, assuming the gauge A^ = and taking 
into account the relation (g|): 

+Nkp{1 - ^) = 0. (24) 



It is easily seen that this equation leads to the field equation (22^ for C -^ +00 and D — Dq — const. 
Authors of pl| prefer to take the continuity equation as the second field equation. By a dimensional 
reasoning the continuity equation (|2^) is generalized to pl| 



;|(p(-n + 4(-)^ = 0. (25) 
at qq at ao 

Now, a qualitative behavior of this toy model is obtained via a first integral of motion. From the 
Lagrangian ([l7|), the Hamiltonian is defined as 

(£> - 1) C2l)2 (5 p „ „^,n /37V2 , , 

H^-^ -—f^3-e^i — )~i7;DS^('^~^)> (26) 



e^( 


ao 


1^- 


dn 




dC 


dt 




dt 



with 



and taking p and p as source [1-3]: 

d, C^ c,^ ^dD-1^2 P 



(27) 



which leads to 



e^( — )^^-^Z)2 _ /^^^2 2C/D] 

dt^ 2kN \q' D3 2 ^ 

= _£.S2^2C/Dd^ (28) 

2 dt ^ ' 



(^^e -^D^)+pDe «o ( ) = 0. (29) 



2kN dt^ D^ ' -^ ' 2 £) 

Taking the derivatives with respect to D instead of t, we arrive at 



-e ooi:)2)+pe "0 (--_) = 0, (30) 



2kN dD^ 1)3 ' ' ^ ^2 D' 



which is equivalent to 



,D—\ _CD_ . „^ /■ , , ,, CD' 

(— -5-e °o D^) + / dD'p{D')e "0 



2kN^ £|3 



x(i-;^)^0. (31) 



This is interpreted now as the equation of motion for D, having a vanishing total energy, in the 
potential 

UiD) := j\D'p{D')e-'^{\ - -^). (32) 

The kinetic energy term is given by 

The matter content is taken in the form of radiation with the following equation of state 

P^^- (34) 

The kinetic energy is always positive for D > \. Therefore the potential energy has to be negative 
for real dimensions. There is a minimum aX D — 2C for the potential (pi) where it is negative. Its 
behavior at large D is given by 

and for D near zero, assuming the pressure remains finite (nonzero) 

U{D)c^-C\iiD. (36) 

Therefore, U tends to infinity when D -^ +oo as well as Z? — > 0. This means that there are two turning 
points for D (where D = 0), one above D = 2C and the other below it. Appropriate assumptions 
about the matter content of the universe would bring the lower turning point to be Z3 > 1. This leads 
to a model universe without any singularity. Even dimension of the model, being a dynamical quantity 
remains finite. The behavior of the spatial dimension with respect to the time near D — Dq — 3 can 
be obtained by expanding ( |3l| ) using D — Dq + e{t). Assuming 

^' « 1, (37) 



they obtain to the lowest order in e 



which has the solution 



Dn 



.2 2kNDo^po 
-''Do'Po.^to-tr, (39) 



2C{Do - 1) 

where to is the time corresponding to Dq. Using this relation, it has then been shown that for e{to) of 
the order of 10~^ the relative change of space dimension is just about two order of magnitude within 
a time range of about 10 times the present accepted age of the universe [ll|-||l3]. That is, as we go 
backward in time there will be almost no change in the space dimension as far as 100 billion years. 
It sh ould be noted, however, that the change of variable from t to D in (pOh is not allowed (see 
Sec. |3.3.2| ). Therefore, the result we just mentioned are questionable. Besides this, the authors of 
[11-13] in the course of the above calculation use the value C = 600, which is not consistence with 
assuming Dpi ~ 3. For this to be the case we must have C ~ 10^, as can be seen from Table I. 

3.2.2 Modification of the original decrumpling model 

Lima et al. p3[ propose an alternative way of writing the field equations. They bring in the lapse 
function and consider the variation of the Lagrangian with respect to it as the second field equa- 
tion. This corresponds to the 00-component of the Einstein field equations (cf. Eq. ( pO| ) for constant 
dimension): 

^ a 2 ^ _2np_ 

m^a' D{D-l) ^ > 



The two field equations (124) and (MG) are consistent witli the continuity equation (Eq). Now, assuming 
a radiation dominated universe (i.e. p = p/D)^ Eqs.(pl Eq, and Eo) leads to the following evolution 
equation for the spatial dimension D 

D^ = -^^ (41) 

C2(i5-1) ' ^ > 

where A = 2kpqN'^ and po is the energy density of the universe corresponding to Z? = Dq. This 
"energy equation" for D corresponds now to (plf). However, there is an inconsistency between them 



(see Sec. 3.3.2). For D < 1 the kinetic term is again negative. Now, the potential energy is always 



negative for 1 < I? < +cx). Thus there is no turning point in this modified model. 

3.3 Critique of the original decrumpling Model and its modification 

As explained before, there is not a unique formulation for the cosmological models with variable 
space dimension, and no principle from which we could derive unique field equations. Therefore, the 
observational consequences are the only way of looking for validity of the models. Despite this non- 
uniqueness and possible observational viability, we would like to mention some theoretical deficiencies 
of the above models. 

3.3.1 Measure of Action 

The main difficulty in formulating a variational calculus for a theory with variable space dimension 
is the dependence of the measure of action on dimension, making the measure a dynamical quantity. 
As we will see in the next section, this dependence seriously influences the field equations. Models 



discussed up to now |11, 13] do not consider this dependence. Taking into account the measure of the 



integral will also change the equation of continuity (see EqJ69). 

3.3.2 Incorrect time evolution equation of the spatial dimension 

We have seen that two equations (|3^) and ( [4l|) leads to quite different results regarding the behavior 
of D. Lima et al. ||l^ are able to use actually one of the field equations directly, i.e. (|40|). Authors of 
Ref. |ll[ , however, use equation (pi|), which is based on a non- valid change of variables in going from 
Eq.(^9|) to (pOl). Therefore, their equation for the dynamics of D, i.e. (pi|), is not valid. 
The point is better understood if we write the Hamiltonian in its correct canonical form. Substituting 
the momentum conjugate to the space dimension, p/j, in (pq), we obtain for the Hamiltonian 






1 

2C2(D "" 

pDj2g2C/I5)^ (42) 



Now, using 7i as a first integral of motion (see Eq.p7|), we obtain 

Since po and D are independent variables, it is obvious that we are not allowed to change the 
differential variable from t to D. Note that Eq.(M3) is just Eq .(p9[ ) where po is substituted in terms 
of D. Indeed, it is easy to show that Eq.(29) or (Esh leads to (I 



4 NEW MODELS FOR A UNIVERSE WITH VARIABLE 
SPACE DIMENSION 

We have mentioned some of the shortcomings of the original model, regarding the field equations and 
its results. Now, as we have mentioned before even the Lagrangian is not unique. Here we try to 

10 



remedy some of the shortcomings and also write down other possible Lagrangians. 
The matter part of the Lagrangian was written in a form such that p and p could be taken as constant 
in the variation of Lagrangian. Now, for the case of variable dimension this procedure may not work. 
Therfore, we prefer to use another procedure to write down the matter Lagrangian which is due to 



Hawking and Ellis |tl7|. This is outlined in section 4.1. In section 4.2, we formulate new actions, and 
Lagrangians, and then derive the corresponding field equations, using the new matter Lagrangian ( p4[ ) 
and implementing the dynamic character of the measure. 

4.1 Hawking— Ellis action of a prefect fluid in constant space dimension 

We want to obtain the energy-momentum tensor of a prefect fluid in constant dimension, via an action 
formalism, in the form 

Tf"" = {p + p)U^'U'' + pg^"" , (44) 

where U is the timelike velocity four vector satisfying 

go.pU'^UP - -1. (45) 

The fluid can be described by a function /i, called the density, and a congruence of timelike flow lines. 
The fluid current vector, defined by j" = ^U"' , is conserved: 

J";a = 0, (46) 

where " ; " denotes covariant derivative. Taking the elastic potential (or internal energy) e as a 
function of /i the action is written as |0| 

Sm = - I d(^+i)xV^ M (1 + e)- (47) 



The calculations may be simplified by noting that the conservation of the current can be expressed as 

r;„ = -L-^(y^r) = o. (48) 

Given the flow lines, the conservation equations determine j" uniquely at each point on a flow line in 
terms of the initial values at some given point on the same flow line. Therefore, {y/—g)j°' is unchanged 
when the metric is varied. But 

/^' = -9-\{V^r){V^f))ga0, (49) 

so 

2p5p={rj^g^''-rr)Sg^,. (50) 

Now, using the definition of energy-momentum tensor ||23| as 

SSm ^\j d(^+i)xy^r^''%,, (51) 

one gets Eq. (H) for the energy momentum-tensor of a perfect fluid in which 

p = Ai(l + e), (52) 

and 

P = A^^(^), (53) 

are the energy density and the pressure respectively. Hence the action for a prefect fluid may be 
expressed as 

Sm = - I S'^+^^x^^gp. (54) 
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Taking a FRW metric given by Eq.(H), it is easy to see that for a comoving observer whose contravariant 
velocity four vector is 

C/" = (iV-i,0,...,0), (55) 

the energy-momentum tensor is given by 

T^"" = diag{pN-^,pa-^, ...,pa-^). (56) 

In the variation of Sm to obtain the energy-momentum tensor we have used the following relations: 

(5/x = (5a, if a -^ a + 5a, (57) 

a 

and 

d> = 0, if N-^N + 5N. (58) 

The matter action (BJ) is now in suitable form to be generalized to variable space dimension. 

4.2 Definition of our model universe and its field equations 

We are now in a position to write down the complete action for matter and gravity in D space 
dimension. Here our treatments are for each kind of D-dimensional topology (closed, flat, and open), 
in contrast to the earlier works in [11-13] which are for the special case of a flat FRW universe. Using 
the relations (|| ^, ^J) we obtain 

<5'o := Sg + Sm 

a^ aN 

Now, trying again to make the action dimensionless, we have to multiply it by the factor Oq °~ (cf. 






Sec. 3.1). Omitting the constant part of this factor, ag °, and integrating over the space part of the 
action (59|), we reach to 



5* := Squ' 



-Do 

ao Zki\ a a 

N^k, 2aN, , , , 

where Vd is the volume of the spacelike sections coming from the space integration due to the homo- 
geneity and the isotropy of the FRW metric. Assuming D to be an integer, Vd is given by 



Vd 






r(i 



if fc = +l, 



fS^Xc-, iffc^O, (61) 



^2ZL^/(Xc), iffc = -l. 



where Xc is a cut-off and f{xc) is a function thereof (see Appendix A for more details). Now let us 
assume D to be a dynamical parameter accepting any non-integer value. So, the volume of spacelike 
sections, Vd, must be defined for fractal or non-integer dimensions. Since the volume of a fractal 
structure having non-integer dimension is beyond the scope of this paper, we take the expression ( ply 
to be valid for the volume of spacelike sections whether the spatial dimension is integer or not. 
Now, the action (pOl) depends on the second time derivative of the dynamical variable a. We can get 
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rid of it through introduction of a total time derivative, as it was done in section |3.l| . We, however, 
intend to generahze our action to the case of variable space dimension. Therefore we have to decide 
which operation being done first: generaliztion to variable space dimension or substitution of the 
second derivative with a total time derivative, as it was done in the case of the original decrumpling 



model in section 3.1 [O. Substituting first the term a with a total time derivative we obtain 



^i '■= TT- / dH— n-iTl-rr) n U7) 



2k J ^ flo^ dt^ aN ' N -g- 

^^)i±)^ - pNVni-ri (62) 

g^ go flQ 



Omitting the total time derivative in the action (p3), we reach at the following Lagrangian 

1 a ^D{D^^ d^ N'^k 
Zk oq jV g a^ 

-pNVoi-f. (63) 

flo 

In the case A: = 0, this is the same Lagrangian as (ll^) except for introducing the dimensional depen- 
dence of the volume and the special form of the matter Lagrangian. As we are going to work in the 
gauge iV = we could derive first the field equations and then insert the gauge condition, or assume 
first that TV = constant and then vary the Lagrangian. But it is easily seen that both procedures 
lead to the same result. Therefore, we will assume from now on, without loss of generality, that 
N — constant. 

Now, we assume the space dimension Z? to be a variable and use the constraint (|^ in order to reduce 
the dynamical variables to TV and a. The equations of motion are then obtained by variation of the 
Lagrangian (psh with respect to N and g. For the variation of p we do need the following relations 



based on Eqs.(Ha, 49) and 



^■"--ivv^^(^^-C)'^")^'' ^''^ 



.2 1 



M^ - - ^.^,.(^).. [(W.(-)V)(iV^.(-)V)]..„ (65) 

and 

,,dlnVD g .dD D 

d^i^-Ri ,r, +ln— )-— + — ](5g. (66) 

dJJ gg gg a 

As mentioned in ( |58| ) variation of p, with respect to N is zero for constant dimension. It is easy to see 
that this fact is still true in the case of variable dimension. 

Now, using (|6^) and the definition of the energy density and the pressure, as given by Eqs.(|52|) and 
( |53| ) the following equations of motion are obtained: 

1 ,a 2 , k '2kp 



and 



(D - 1 - + - ' + -^)i-T^—nr- 

a a g^ 2C dD 

D(2D-1) D^ , , ., dlnVoD 



2C{D-1) 2Do' ' ' dD C 
a 
G go 



^. In — + 1)] = 0. (68) 
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Using Eqs.(p^,pq), one easily gets the continuity equation for our model with variable space dimension: 

;|(p(-)^^z5)+P;|((-)^^z.) = 0. (69) 

at flo at gq 

This continuity equation can be integrated for the case of dust or radiation to obtain the energy 
density as a function of time. For radiation era, p = -^j we have 

p^Poe '°o ^( )Da e -""o , (70) 

and for matter era, p = 0: 

p^pae ^oo '-—. (71) 

VD 

Using Eqs.(ra, pTj ), the evolution equation for the spatial dimension is derived as 

However, if we assume the time variability of D first, substitution of the second time derivative by a 
total time derivative leads to the following action 

1 f , (r, '^ ,DVd a . a ]j 2DaVD , a ,^ 



Sii ■■= - dt{2-{—^-{-)^) -ii(-) 

Ik J at I\ a ag al\ oq 

DVd at) a a jj 2DDa a^dVo 



N a ao ao N a ao dD 

N a a'^ ao 

-pNVD{-f}. (73) 

ao 

Neglecting now the total time derivative, we reach at the following Lagrangian 

VD,a jj2ba 2Dab a D{D-l) 
" '^ ~^^a^^ ^~^^ aN '^^^ TV 

,,a,n N'^k. 2Db a , a .T-idlnVo^ 

x((-) —) + -^r ( — ) ,n } 

a a^ Al a ao dU 

-pVdN{-)''. (74) 

ao 

This Lagrangian is more general than L/. Putting Z? = in L// we obtain the Lagrangian Li which 
can still be considered as a Lagrangian with variable space dimension. 

Using again the dimensional constraint (j^) and Eqs.(p2| , [53|) and (pq), variation of Lu with respect 
to N and a leads to the following equations of motion 



-, . r, kD(D-l) 

^„> r^r2D2 T^ 1 2D 2D2rflnVm' \"^ I 



and 



N^'a' n[^-D-l---^---^^-^ 



h,2D 2D^ ^ a 2D^ dlnVo 
a C C ao C dD 

>,,, £|4 a,, , 4D3 3D^ 51)2 ^ 

+ (-)'h7^(ln— )' + i--^ + ^5?^ - ^^1^ — 
a C^ ao C^ zC 2o ao 
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4 



a 



dlnVoAD^ 3D^ 5L»2 2D „ 
dD ^~CP 2C '^ '2C ^ "C2" ^ ~q' 
D^ d'^Vo 31)2 2i:>2 5£, (£»-!)(£) -2) 



C2Vb rfL'^ ' c C2 2C 2 ' 

fciV2 n(2L>-l) L>2(D- l)dlnFz3 
~^^ 2C ^ 2C^ dD 

(£»-!)(£) -2) D^{D~l)^ a, ..2 /^, « 

-■^^ ^^^ ^ H ^-7; In — ] + N^Kp{— In — 

2 2(7 ao C ao 

DdlnVo 



In the limit C — > +00 and D = Dq = const, Eqs.(]6^, pSl Uq, 76) reduce to the standard Friedmann 
equations ( p0| , p2| ). Using Eqs.dTS, ^q), one easily gets the continuity equation for our model with 
variable space dimension. Now, using Eqs.(1q, U^), the evolution equation for the spatial dimension is 
derived as 



^2 _ N^D^[2Kp-kD{D-l)5- 

^2r 2-D^ jj ^ 2D 2D^ dlnVp ] 



-^ ^Tr2J)^ n 1 2D 2D^ dlnVp T- ^' ' > 

D„ ^ ^ C C dD 



From Eqs.(|72|) and (|77|), it is easy to see that in the limit C -^ +00, we have D ~ 0, which corresponds 
to the standard universe model with constant dimension. To understand the dynamical behavior of 
dimension in the models discussed we have to substitute p and Vd from (|7(]|-|7l|) and (^) in (JT^) or 
( [771 ) respectively. In the case of fc = it is, however, seen easily that there is no turning point in 
dimension and both models can be contracted up to a = S corresponding to £* — + +00. The discussion 
can be simplified for a de Sitter-like model which we will discuss in the next section. 

4.3 de Sitter— like model 

Consider a universe dominated by the constant vacuum energy, p\, corresponding to a cosmological 
constant A and the equation of state pA ~ —pK- Note that these relations are satisfied by the continuity 
equation (pOf). Substituting 

in ( [72[ ) and ( [77] ) , we obtain evolution equation of the spatial dimension for the Lagrangian Li 

_i\r2jj3 
D' = c'HD-l) ^^^^ ' l)S-'e-^^/^ 2A] ^ -Vj{D), (79) 



and for the Lagrangian Lji 



-.2 



N^D^[D{D - l),5-2e-2C/i3 _ 2A] 



^" = ~^2r2D^ 2D^ dlnVn 2D ^ ^1 = "^^K^)' (80) 

L Dq C dD C J 



where we have set fc = +1. As was mentioned in Sec. 3.2 , these equations may be interpreted as "energy 
equation" for D, having vanishing total energy with the potentials V/(D) and Vii{D) respectively. 
Both models have a "classical" turning point at Dt, where D vanishes. This is the point of maximum 
dimension of space corresponding to the minimum radius of the universe ut (see Fig. 5). We obtain 
this dimension as the solution of the equation 

£>(£>- l)5-2e-2C/i5 _ 2A = 0. (81) 

In what follows we will choose for the cosmological constant 

A = 3G"^ =3/pr^. (82) 

Eq.(pl|) has to be solved numerically. Table III shows some interesting Dt values. There is a minimum 
dimension for both models which is reached asymptotically. This minimum dimension for model I is 
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D = 1 corresponding to a maximum radius a/„i = 6e'~' . To study the asymptotic behavior of Vii{D), 
we need to calculate the term '^ ^dX" ■ From (bH), for fc = +1, we obtain 



dinVp 1 



D + 1 



)), 



(83) 



where ■0 is the logarithmic derivative of the gamma function, ip — y?V)' called Euler's Psi function. 
Substituting ( |8^ ) into (|8(]|), it is easy to see that the potential V//(D) has an asymptotic behavior at 
Diirm corresponding to a//„i, being the solution of the equation 



20"^ 1)2 n + 

— - In TT — w 

Do C ^ ^^2 



i,,-^-.-i^o. 



(84) 



Some numerical values for Dum and a//„j are also given in Table III. A generic example of the 
potential Vii{D) is plotted in Fig. 5. It exhibits the classically allowed region {D^ > 0). Therefore 
in a de Sitter-like universe, taking the model II, for the dimension of space we have Djjm < D < 
Dt- Dt, which corresponds to the minimum contraction in the de Sitter-like model, can have very 
different dimensions, as given in Table III. There we can also see that the minimum scale factor, 
ar, corresponding to Dt can be up to about 10 times the Planck length, assuming Dq = 3. The 
assumption Dq = 2 leads to almost the same result with values for Dt a bit smaller than those of 
Table III. The maximum radius aum, corresponding to the minimum dimension Dum which is about 
2 as can be seen from Table III, is many order of magnitudes bigger than the present radius of the 
universe. 

The behavior of Vi{D) is qualitatively the same as V7/(D), except that the minimum dimension is at 
Djim = 1. In the next sections on quantum cosmology of our de Sitter-like models we will come back 
to the details of this behavior. The detailed behavior of the potentials and the field equations leads 
to interesting cosmological consequences which we will discuss in a forthcoming paper p3. 
TABLE III. aT/lpi,DT,log{aiijn/lpi),Djjjn,log{airn/lpi) for A: = +1 and for interesting values of C 
corresponding to Dpi — 3.001, 4, 10, 25 and Dq — 3 (cf Table I for corresponding values of S). 



c 


ar/lpi 


Dt 


log(a//„i/;p;) 


Dllrn 


log(a/,„//p;) 


1.2588 X 10*^ 


1.0 


3.001 


7.8160 X 104 


2.0 


3.6452 X 10^ 


1.6788 X 10^ 


1.4 


3.996 


1.5079 X 10^ 


2.189 


5.4682 X 10^ 


5.9957 X 10^ 


3.0 


9.782 


9.2644 X 10 


2.194 


2.3435 X 10^ 


4.7693 X 10^ 


8.0 


22.42 


8.6034 X 10 


2.196 


1.9884 X 10^ 



5 WAVE FUNCTION OF OUR MODEL UNIVERSE 

We are now interested in quantum cosmological behavior of our model universe and its differences to 
the Wheeler-DeWitt equation and its solutions for the de Sitter minisuperspace in 3-space. We will 
therefore use the canonical approach of quantization due to DeWitt and Wheeler ||2|, ^. To do so, 
let us first briefly review the canonical quantization of the universe with constant spatial dimension. 

5.1 Brief review of the tunneling and the Hartle-Hawking wave function 

According to the quantum approach to the standard cosmology p^, Oq, p9 , a small closed universe 
can spontaneously nucleate out of "nothing" , where by "nothing" we mean a state with no classical 
space-time. The cosmological wave function can be used to calculate the probability distribution for 
the initial configuration of the nucleating universe. Once the universe nucleated, it is expected to go 
through a period of inflation. 

Let us now illustrate how the nucleation of the universe can be described in the simplest de Sitter 
minisuperspace model. In this model, the universe is assumed to be homogeneous, isotropic, closed, 
and filled with a vacuum of constant energy density pA- The universe should be closed, since otherwise 
its volume will be infinite and the nucleation probability would be zero. The radius of the universe a 
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is the only dynamical variable of the model, and the wave function ^(a) satisfies the Wheeler-DeWitt 
(WDW) equation: 

H^^ia) = 0, (85) 

where H is the corresponding Haniiltonian of the model. In our de Sitter case it is written as 

d2 9^2^2 



Here 



&-^(1-^a'«')]*(«)=0- (86) 

2 _ SnGpA _ A 



The value A = 3G~^ corresponds to l\ — Ipi. We have disregarded the ambiguity in the ordering 
of non -commuting operators a and d/da, which is unimportant in the semiclassical approximation 
we are using here. Eq. ( pq ) has the form of a one-dimensional Schrodinger equation for a "particle" 
described by a coordinate a(t) having zero energy and moving in the potential 

^(«) = ^(1-?aV). (88) 

The classically allowed region is a > /a, and the WKB solution of Eq.(p6|) in this region is 

*±(a) = [p(a)]"^/^exp[±i / p{a')da' T iTr/4], (89) 

Jia 

where p(a) = [— C/(a)]^/^. The under-barrier (a < Za) solutions are 

^±(«) = b(«)r'/' exp[± / ' \p{a')\da']. (90) 

J a 

The classical momentum conjugate to a is pa = ~ 2GN - ^^^ a^ Ia ^e have 

Pa*±(a)«±p(a)*±(a), (91) 

where pa = ^id/da. Thus \l/_(a) and ^+(a) describe an expanding and a contracting universe, 
respectively. According to Vilenkin, the tunneling boundary condition BQ] requires that only the 
expanding component should be present at large a. Therefore we obtain for the Vilenkin tunneling 
wave function ^y 

*y(a>^A) = *-(a). (92) 

The under-barrier wave function is then found from WKB connection formula: 

*y(a<^A) = *+(fl)-^*-(a). (93) 

The growing exponential ^_(a) and the decreasing exponential ^+(a) have comparable amplitudes 
at the nucleation point a = Ia, but away from that point the decreasing exponential dominates (see 
Fig. 6). Therefore, the nucleation probability can be approximated as [|8| 

^.,,*v'(?A),2 ,_r , r'^...'^u.n _„,... Stt 



exp[-2^'b(a')Ma']=exp(-|^). (94) 



' *i^(0) ' 
The Hartle-Hawking wave function satisfying the no-boundary boundary condition takes the from 

ii 

*//ff(a<^A)-*-(a), (95) 

for the under-barrier wave function and 

*//ff(a>/A) = *+(a) + *-(a), (96) 

17 



in the classically allowed range. This wave function describes a contracting and re-expanding universe; 
under the barrier '^Hnia) is exponentially suppressed (see Fig. 6). 

Using the anti-Wick rotation for Euclideanization of a Lorentzian path integral, Linde |S| suggested 
his tunneling wave function. Linde's wave function has only the decreasing exponential ^+(a) in the 
non-classically allowed region: 

^L{a < Ia) = ^+{a). (97) 

The continuation to the classically allowed region gives 

*L(a>ZA) = ^[*+(a) + *-(«)]. (98) 

Different probability densities are assumed to correspond to different boundary conditions implied by 
the Hartle-Hawking "no-boundary" proposal [2^ , the "tunneling" proposal of Vilenkin [^ ^ , and 
that of Linde p9| |. In particular, the nucleation probability for instanton-dominated transitions is 
assumed to be 

' ' \exp(2S'B) for fL, *v, ^ ' 

where Se is the Euclidean action of the instanton. Problems concerning the interpratation of the 
wave function of the universe are studied in many articles |^3[ . For a recent discussion of problems 
associated with defining the initial cosmological wave function see e.g. Ref.[|4). The debate about 
the form of the wave function of the universe has recently intensified by Vilenkin ||35| , Hawking and 
Turok 1^, and Linde 0; see also Refs.|3|, |3|. 

5.2 Wheeler— DeWitt equation in the universe with dynamical spatial di- 
mension 

Turning now to the canonical quantization of our model with variable space dimension, we consider 
the simplest case of the de Sitter minisuperspace. We therefore take k = +1 and p = pa — g^j in 
the Lagrangian Lj and L//. The corresponding Hamiltonians to Lj and Lu can be written in the 
form 

Hi ■- NHi 

:= piao, - Li 



na^Pia^ _ VdDJD-I) a. 



VdA, a ^ 



i-r}, (100) 



and 



Hi I := NHii := piiaO- - Ln 

Ka^PlIa 



2 



N{- 



ADV (^\D\Dl _ D^ dlnVp _ D _ D+1 
^^^D{ao> I- Do C dD C 2 



- ^^^'r" (-)" + ^(-n^ (101) 

where Vd is the volume of D-sphere S^ (cf. Appendix A and Eq|6l|). The canonical momenta pia 
and piia are defined as 



Pla 



dLi 
da 

Kl\ a^ ao 
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and 

dLii dLii dD 



Plla 



da db da 

VoDd a rj 2D^ 2D^ dlnVo 
~ KNa^ ^'^' ^Id'o C dD 

^-D-1]. (103) 



The above systems can be quantized by the assignment 

P/a^ma^--ao-'^"^|-a'^|-, (104) 

ai da da 

where q is the factor ordering parameter. Note that we have introduced the coefficient ao^^^° in ( |104[ ) 
for the kinetic term in WDW equation to have the right dimension. The WDW equation is obtained 
by applying the classical constraint 

^^' = 0, (105) 



5N 
SHii 



SN 
on the wave function ^/(a) and '^jj{a) respectively: 



= 0, (106) 



^'^'^"^ = ^^d'a^'d'a ^^ (^^ 

X (^^i£-^-2A)}^iia)^0, (107) 



and 



D{D-l) 2D^ 2D^ dlnVD 

"" ^ a^ ^^>^^~~C dD~ 

2D 
- — -i?-l]}*//(a)=0. (108) 

In the semiclassical approximation the ambiguity in the ordering of a and d/da can be ignored. 
Therefore, we will assume q ~ 0. The corresponding potentials for ^/(a) and ^//(a) are 

. (2^-2A,. (109) 

2D^ dlnVo 2D 



C dD C 



D-l]. (110) 



It is worth noticing that in these equations 13 is a function of the scale factor, a, according to the 
constraint (ph. It can be easily shown 

Ui{a) -^ U{a), as D = Dq = 3, 
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and 

Un{a) -> U{a), as D ^ Dq = 3, C ^ oo. 

Note that C appears explicitly in Uii[a). This is due to appearance of the time derivative of D 
in the Lagrangian L//. The shape of the potential Ui{a) and Unia) determine how we impose the 
appropriate boundary condition and evaluate the wave function of our model. Using ( |6l[ ) for fc = +1 
and ( |83| ) and the dimensional constraint (^, the potentials Ui{a) and Uii{a) can be written in terms 
of D only: 



Ui{D) 

and 

Uii{D) 



(£>(£» -l)(5-2e-2C/i5_2A), (HI) 



4ao'^°7r(^+i)D<5-2e'^(i-*-*) 



^Hn^)) 



2 

2D^ 



(£>(£» -l)5-2e-2C/s_2A)( 



Do 
- ^(ln--V^(^))-^-i?-l). (112) 

Let us now obtain the zero points of Ui{D) and Uii{D). Using the following asymptotic expression 
for r and '0 function as D ^ +00 

r(i?)-e-^i?^-*^/2^(l + ^ + ...), 
and 



e-^i?^- 


"^(^+12Z.+ 


^InD- 


1 V^ B2k 

2D ^ 2kD^'' 

fe=0 



where -82*; are the Bernoulli numbers, we obtain for D -^ +cx) 

U,[D) ^ ^^(2^)^i^(4-z,)gD(i-|§) ^ 0^ 

and 

Therefore, the potentials Ui{D) and Uii{D) both tend to zero as _D ^ +00. According to Eqs.(^-p4) 



and using Eqs.(lll 112) three other zero points of Ui{D) and Uii{D) are obtained. Those of Ui{D) 
are 



D- 


> +00 


for a - 


-5, 


Dt 




for a ~ 


= AT, 


1 




for a - 


= a/m, 







for a - 


-> +00 



^ = <^ -T -. " - -T, (^^3) 



For Uii{D) we have the same zeros except the third one which is at Z? = !)//,„ corresponding to 
a = ajirm whose corresponding values are given in Table III. 

We are now in a position to compare these potentials with the potential of the standard de Sitter 
minisuperspace depicted in Fig. 6. To make the comparison as simple as possible, we draw the 
dependence of our potential Ui{a) against the scale factor a. As we see from Fig. 7, the behavior of 
the potential [// (a) for small scale factors of the order of the Planck length is similar to the potential 
U{a). Instead, for large scale factors we see a completely different behavior. Here we have two 
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potential barriers with different heights depending on the value of C. For example, for C = 1678, 

corresponding to Dpi — 4, the ratio of their heights is of the order of 10^''®. 

The height of the potential barrier in the region (5 < a < ot is also interesting to be compared to that 



of the potential U{a). From Eq.(88), it is easy to see that the height of the potential barrier of U{a), 
in the region < a < Za, is of the order G~^, or the square of the Planck energy, which is of the order 
of lO'^^Gey^. In contrast, the height of the potential barrier of Ui{a) in the region 5 < a < ar ^s 
about 10-i58Gey2 foj. c = 1678 and IQ-^'^^^GeV'^ for C = 477. Therefore as C increases, the height 
of the barrier increases so that for C ~ 10^ it is of the order G~^ comparable to the potential barrier 
U{a). The behavior of Uii{a) is similar to Ui{a). 

5.3 Wave function of the universe with dynamical space dimension 

Here, we are interested in the solutions of the WDW equation in our model universe. Let us first use 



the semiclassical approximation to solve Eqs. (107) and then Eq.(108). The classically allowed region 
is ar < a < aim, and the WKB solutions are 

*/±(a) = [pM]'^'^ exp[±i / pi{a')da' T J^/4], (114) 

J aT 

where pi{a) = [— C//(a)]^/^. Under the potential barrier in the region 6 < a < ar the WKB solutions 
are 

*7±(«) = b/(a)r'/'exp[± / " \pi{a')\da']. (115) 

J a 

Finally, the solutions for the other under-barrier region a > aj^n cire 

*/±(a) = b/(a)r'/'exp[± /" \pi{a')\da']. (116) 



From Eq.( |l02D we see that the classical momentum conjugate to a, pja, is proportional to —a{D — 1). 
For gt ^ a < ajm, Eq.(|9l|) can be used for ^/±: 

Pa*/±(a) ~ ±p/(a)*/±(a), (117) 

Thus ^/_(a) and \l//_|_(a) describe an expanding and a contracting universe, respectively. Now, we 
have to impose boundary condition in order to specify the wave function of our model uniquely. 
We choose, in accordance with the tunneling and no-boundary boundary condition, the following 
boundary condition as the scale factor tends to infinity: 

lim */(a) -> 0. (118) 

a — *+oo 

Another boundary condition must be imposed at the turning point a = 6 where the semiclassical 
approximation is bound to be wrong. Taking into account the normalization of the wave function 
pi| , we can write it as 

*/((5) = 1. (119) 

This satisfies the regularity of the wave function, i.e.| ^| < oo, as in the case of Hartle-Hawking p8| 



and Vilenkin |3l|] boundary conditions. According to (lli) the growing exponential ^/4-(a) should be 
absent under the barrier in the region a > aim '■ 

*/(a>a/m) = ^7-(a). (120) 

In the classically allowed region, the wave fvmction is found by using the WKB connection formula at 
the turning point a = aim |4^ : 

^i{aT < a < aim) 

= 2[pi{a)]-'/^ cos{ r^ pi{a')da' ~^). (121) 
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This may also be written as 



"^iiiiaT <a< aim) = 2[p/(a)] ^/^ 



X sin( / pi{a')da' — / pi{a')da' + —) 
= 2[pi{a)]-'/^sm{ T"^ pi{a')da') 



X cos( / p/(a')(ia' — — ) — cos( / pi{a')da') 



sin(/ pj{a')da' - -)}. 



(122) 



Using (114), this gives 



^/(oT <a< aim) = 'Si-{a)e ■''■■t 



+ */+(a)e 



-i{j°'''" pi[a')da' -^) 



(123) 



Note that the dynamical character of the space dimension makes the wave function in the classical 
region a mixture of expanding and contracting part. It is in contrast to the standard case of con- 
stant space dimension where we can choose the wave function to be just expanding as assumed in the 
Vilenkin's tunneling boundary condition pQ]. The presence of the expanding and contracting compo- 
nent in the classical region is due to the special form of our potential (see Fig. 7). The under-barrier 
wave function is found by using the WKB connection formula at the turning point a — aT'. 

/■aim 

^i{5 <a<aT)^ \pi{a)\-^/^{sm{ / pi{a')da') 



X exp(- 



\pi{a')\da) -f cos( / pi{a)da 



X exp( / \pi{a')\da')}. 



(124) 



From ( 115 ), this may be written as 

raim 

■^i{S <a < ar) = sin( / pi{a')da')^i^{a) 

^ aT 
raim 

+ cos( / pi{a)da)^i+{a). 



(125) 



The wave function ^/(a) is schematically represented in Fig. 7. Similarly, we obtain the following 
wave functions for the potential C///(a): 



(126) 



*//((5 <a<aT) = sin( / pii{a')da)'^ ii-{a) 

J aT 
raiim 

+ cos( / pii{a')da')^ii+{a), 



^ii[aT < a < aiim) =^ii-[a)e •''■^T 

-i(J °-'^™ pii{a')da'-^) 



+ *//+(a)e 



and 



*//(a > aiim) = *//-(«). 



(127) 
(128) 



22 



The above wave functions are different from those of Vilcnkin, Hartle-Hawking, or Linde. But, as 



we will see, there are similarities between our general wave functions given by (12C, 123, 125) for Ui 



or (126, 127, 128) for Un and those of Vilenkin, Hartle-Hawking, and Linde. We notice first that 
our wave functions are real for every dimension, even for the limiting case of C — > +oo or 13 ^ Dq. 
Therefore, we conclude that our general wave functions are not of Vilenkin's type which is complex. In 
the classically allowed region, '^hh and ^l are both a superposition of the expanding and contracting 
terms with equal coefficients. In contrast, in our general wave functions the coefhcients of expanding 
and contracting components are different. Looking at our general wave functions, wee see that both 
decreasing and increasing exponential terms are present in the non-classical range oi 6 < a < ax- 
These means that in general, in this region, we ex pect our wave functions are a superposition of '^hh 
and ^L. Requiring the decreasing term, \E'/+, in (|l25|) to vanish we are led to: 



aim 1 

p/(a')da' = (2n+-)7r, (129) 

where n is an integer number. This requirement makes the coefficients of the expanding and contract- 
ing part of our wave functions in the classical region to be equal. Hence, assuming the relation ( |129[ ) 
our general wave functions behave as "^ hh in the region S < a < ajm- 

Now, if we require the increasing term "^j- to vanish under the potential barrier in the region of 
(5 < a < Or, we obtain 

pi{a')da' = 2rmT, (130) 

/ aT 

where m is an integer number. This leads to the following wave functions: 

*^'"^((5<a<aT) = */+(a), (131) 

^'^;^\aT <a< airn) = i [*/+(«) - */-(«)], (132) 

*[™^(a>a/,„)EE^z-(a), (133) 

which are similar to the Linde's wave function. ^}' behaves as ^l in the range of 6 < a < ar- But 

in the classical region, i.e. ax < a < ajm, there are differences to ^l- We call ^^ the modified 

Linde wave function. 

Now, we are interested in the form of our wave function in the limiting case of constant space dimension, 

i.e. C — > +00. In Appendix B, we show that as C ^ +oo the relevant integral behaves in the following 

way: 

raim 3 2C 

hm / pjia')da'= lun -JL_ ^ +^. (134) 



Substituting this limiting behavior in (123) and (125), we see that the relevant terms are not well 



defined. But we may assume cases in which C accepts those values corresponding to the relation 
( |129| ) or ( |13C| ). It can easily be seen that in this limit \&/_(a), \E'/+(a) and ^/_(a) tend to ^-(a), 
\I'+(a), ^^(a), respectively. Now, taking the limit of C to infinity corresponding to n ^ +oo, our 
wave functions approach to the Hartle-Hawking one: 

lim ^["^((5<a<aT) = *HH(0<a<ZA), (135) 

lim *["'(aT <a<a/m) = *HH(a>^A). (136) 

n — >+oo 



If the limiting behavior of C in the relation (130) corresponds to to — *■ +00, then we obtain the 
modified Linde wave function: 



,,(™)^ 
,(")^ 



lim *^"^((5<a<aT) = *L(0<a<ZA), (137) 



lim *;^"^(aT < a <a/„0 =«[*+(«)- *-(a)]. (138) 
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It should be noted that in this hmiting case, the second barrier in the large scale factors is removed 
and our potential Ui{a) tends to the standard potential U{a). Similarly, we can show that the same 
is true for the wave function '^ij{a). 



We conclude that assuming the boundary condition (118) and using WKB approach we always have 
contracting and expanding terms in the classicl region, i.e. Vilenkin wave function can not be derived 
from our general wave function. This is due to the form of our potential which has a second barrier 
at large scale factors, independent of the C-value. Therefore, we may say that the Vilenkin wave 
function is structurally unstable with respect to the variation of dimension. 

5.4 The probability density 

It is interesting to find the probability density in our model universe and compare it to its value of 
the de Sitter minisuperspace in 3-space. From Eq.(M), we can calculate the probability density 

Vi^exp{~2 \pj{a')\da'), (139) 

Js 

for the potential barrier of Ui{a) in the region 6 < a < ar- As before, we take A — 3G^^. The integral 
can be calculated numerically for different values of C: 

rexp(-2.72) if C= 1258800, 

^ "" \ exp(-2.28 X lO-'^O) if C = 1678.8. ^ ' 

The probability density decreases as C increases, i.e. as the height of the potential barrier increases. 
Consequently, in the limit C —f +oo, corresponding to the constant space dimension, the probability 
density has its minimum value given by Eq.(|9J): 

lim 7^/ =P-exp(-— -) =exp(-7r). (141) 

c—>+oc GA 

This is the same as the tunnelling probability proposed by Vilenkin and Linde. There is recently a 
controversy between Hartle, Hawking, and Turok on one side and Vilenkin and Linde on the other side 
on the sign of the action in the probabity density (Pj|-[p9|, pq|). Our calculation of the probability 
leads to the proposal of Vilenkin and Linde. There is another barrier in the region aim < a < +oo 
for the potential Ui, as seen by Fig. 7. The probability density for this potential barrier: 

r + co 

'Pi^expi-2 \piia')\da'), (142) 

is less than the probability density of the potential barrier in small scale factor. This is due to the 
height of the barrier for large scale factor is much more than the height of the barrier in small scale 
factor, see Fig. 7. The probability density 

Vh ~ exp(-2 / ^ \pii{a')\da'), (143) 

Js 

for the potential barrier of Uu (a) in the region 6 < a < ax and 

A/~exp(-2/ \pii{a')\da'), (144) 



for the potential barrier in the region a//„i < a < +oo can similarly be calculated and leads to similar 
results like Ui{a). The existence of this barrier for large scale factors and the behavior of the potential 
in that region has cosmological consequences which will be discussed in a forthcoming paper l24| . 
Assuming A = 3G~^ ~ lO^^Gey^, the Euclidean action Se = (— ^) takes the value (— tt). This large 
value for A-term is estimated in modern theories of elementary particle. Astronomical observations 
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indicates that the cosmological constant is much less than this vahie. Assuming A < IQ^^^GeV^, the 
Euchdean action is very large and negative {Se ~ —10^^° typically). This gives the probability density 
proportional to P oc exp(— 10^^*^). Therefore, the probability density in the de Sitter minisuperspace 
depends on the assumed value for the cosmological constant. From the above discussion we conclude 
that for a given value of the cosmological constant, the probability density in our model is much more 
than the probability density of the de Sitter minisuperspace in 3-space. Therefore, the creation of a 
closed Friedman universe with variable spatial dimension is much more probable than the creation of 
a closed de Sitter universe in 3-space. 

So in our model, the universe tunnels from the first turning point, a — d to the second turning 
point, a = ax, expands up to the third turning point a = ajm (or aum for Un). It may then starts 
tunneling through the potential barrier in the region a > aim (or aum)- Seeing the process classi- 
cally, the point a/„i (or ajjm) may act as a turning point at which the universe begins to contract |}42| . 



6 CONCLUDING REMARKS 

The idea of the space dimension to be other than three is relatively old, at least since the work of 
Kaluza and Klein. But it seems that the variability of space dimension is a relative new idea. We 
accept it as a viable alternative and formulate a cosmological toy model based on a Lagrangian to look 
for its consequences. In our formulation, there is a constraint and a scalar parameter, C, accepting 
every positive value. The limiting case of C — > +00 corresponds to the standard model of constant 
space dimension. Taking the present spatial dimension, Dq, to be equal to 3, and assuming the space 
dimension at the Planck length. Dpi, to be about 3, or one of the values of 4, 10, or 25 coming from 
the superstring theories, the corresponding values of C are of the order 10^, 10^, 10^, 10^, respectively. 
For Do = 2 as a fractal dimension for matter distribution in the universe, coming from cosmological 
considerations, we obtain different values of C corresponding to Dpi = 3, 4, 10 or 25. 
After a critical review of the previous works on the decrumpling universe model, we generalized the 
Hawking-EUis action for the prefect fluid universe in 3-space to the case of variable space dimension. 
Using this generalized action for a perfect fluid, we have then formulated the action for a FRW 
universe filled by a prefect fluid and allowing the space dimension to be a dynamical scalar variable. 
We have seen that this generalization is not unique. In contrast to the earlier works, we have taken 
into account the dependence of the measure of the integral on the spatial dimension. The Lagrangian 
and the equations of motion and also the time evolution equation of the spatial dimension are then 
obtained. 

Using the time evolution equation of the spatial dimension, we have obtained the classical turning 
points of our model universe. It turns out that the space dimension could have been up to infinity at 
the beginning of the expansion phase of the universe, where the universe has a minimum size equal 
to S. This is the only singularity we encounter in our model. 

We then turn to quantum formulation of our toy model and write down the Wheeler~De Witt equation 
for our universe. Imposing the appropriate boundary condition at a — > +00 and using the semicalssical 
approximation we have obtained the wave function of our model, which is real but more general than 
its standard counterparts. As a result, the Hartle-Hawking wave function is obtained for special values 
of the constant C. Other values of C lead to a modified Linde wave function. In the standard limit of 
constant space dimension we may have a general, Hartle-Hawking, or a modified Linde wave function, 
but that of Vilenkin is not obtained. This is due to the special shape of the potential which has 
always a second barrier for large scale factors and only vanishes in the limiting case of C — > +00. We 
have therefore concluded that the Vilenkin wave function is, in respect of varing the space dimension, 
structurally unstable. 

Finally we have calculated the probability density for the potential barrier in small scale factor. In 
the limit of constant spatial dimension, this probability density approaches to its corresponding value 
of the de Sitter minisuperspace as suggested by Vilenkin and Linde. In general, the probabilty density 
for our toy model is larger than that of the standard quantum cosmology. 
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Our toy model shows that it is possible to formulate a Lagrangian for a cosmological model with 
dynamical space dimension. It removes the usual singularities in physical quantities like scale factor 
of the universe because the minimum scale factor of the universe is 6. It gives us a model of how a 
probable space dimension of 10 or 25 at the Planck length can be incorporated in the usual picture 
of (1 + 3)-dimensional space time and what its consequences are. We are currently studying other 
cosmological consequences of the model, specially those which can be compared to obsevational data. 
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APPENDIX A: Volume of spacelike sections 

The D-dimensional space metric is defined as (see Eq. 0) 

+ sm^ ei...sm^ eD-2d9%_i], (145) 

where 

O^e.^TT i = l,...,D-2 0<6lD_i<27r. 

The radius F{x) is expressed as 

{sinx ii k = +1, < X < ti", 
X iik = 0, 0<x<oo, (146) 

sinh X ii k = —1, < X < oo. 

We know consider the different cases fc = 0, ±1, separately: 

i) fc = +l: 

Vd is the volume of a unit D-sphere, S^: 

Vd : 



/^in^- 

"'0 


'xrfx / dOD-i 
Jo 


D-1.D>2 

n 


[ d0o-mSin™-i 


m=2 


h 


2,r(^) 





'D-m 



ii) k = 0: 

In this case Vd is infinite. Introducing a cut-off, Xc, as a very large number, Vd has the form 



(147) 



Vd = 


/•Xc /-Ztt 

/ X^'-^dx / dOD-i 
Jo Jo 




D-1M>2 ^ 


X 


\\ / d0i3_,„sin'"-i^ 




m=2 -^0 




^^D^D/2 




r(# + i)' 


iii) fc = -1: 




Again, using the case of xc, we obtain 


Vd = 


/■Xc r2iT 

/ sinh^-^xrfx / dOD-i 
Jo Jo 



D-m 



(148) 
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D-1,D>2 



„ n iJ 



X 11/ dt/D-m sm™ 1 6li3_^. 



m=2 



/(Xc). (149) 



r(f) 

where 

f(-^\^) 2--^"'-i A^fe=OV -y V fe ; 2m-2fe ' "" --"■■" ' -' (^r.o.\ 

jy?ic)—\ 1 ,^„j_i T \k(2m-l\ cosh(2m-l-2fc)Xc-l ■ -• t^ r, l-^^^j 



1 V^'" / 1 \k(2rn\ sinh(2m-2fc)xc if H — 9Tn ^ 1 

25^^T^z^fe=ol~-^J UJ — 2;;r:2fe — ' iti^-zim + i, 

1 Y^m-1/ T Nifc /'2m-l\ cosh(2m-l-2fc)Xc-l i-P n — O:,- 

W^^r^ 2^k=oy~^) [ k ) 2m-i-2fc ' ^i ^ - ^^ 



where m is an integer number. For details of the about integration see p6| . In general xc tends to 
infinity so that for k = 0,-1, the volume of the spacelike sections are infinite. 

It should be emphasized that the above expression for Vu is valid for integer dimension. Volume of 
fractal structures having a non-integer dimension is beyond the scope of this paper. 

APPENDIX B: Limiting behavior of our wave function for C -^ +00 



Our goal is to obtain the behavior of integrals in the wave function ( |l23| ) and (125) as C 

/•aim 

lim / pr(a)da. 
We calculate it by expanding pi(a) in inverse powers of C: 



, , Stto ,0? , f(a) ^^, 1 , /,^,x 

where we take Dq — 3 and 

where ■0(2) — 0.42. From Table III, we know that for C ^ 10^ or so, the value of a^ is of the order of 
the Planck length. Therefore, we can also expand ay in inverse powers of C: 

aT = VG+^+0(ij) + ..., (153) 



where A is a constant. From Eqs.(151)-(153) we obtain 



^^"'" pj{a)da = \{-l+ '^f" + ^ + 0{^) + ..., (154) 

where Ai is a constant. As previously mentioned, the relation between ajm and C is 

aim = aoe~. (155) 



Here ag is the scale factor corresponding to Z? = Dq. Inserting (|155|) into (154), we are led to 
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FIGURE CAPTIONS 

FIG. 1. The values of C as a function of Dpi when 4 < Dpi < 26 and Dq taking the values 2.9 
(dashed line), 3.0 (solid line), 3.1 (dotted line). Using Eq.(||) it is easy to show that if we take Dq — 3 
for Dpi — 3 and Dpi -^ +oo, we have C — > +oo and C — Do In ^ — ~ 419.7, respectively. 

FIG. 2. The values of C as a function of log Dp; when 3.6 < log Dp; < 6.2 and Dq taking the 
values 2.9999 (dashed hue), 3.0000 (sohd line), 3.0001 (dotted line). 

FIG. 3. The values of log | log j^| as a function of Dpi when Dpi ~ 3 and Dq taking the values 
2.990 (dashed line), 2.995 (solid line), 3.000 (dotted line). 

FIG. 4. The values of log( j^) as a function of Dpi when 4 < Dpi < 26 and Dq taking the values 
3.1 (dashed line), 3.0 (solid line), 2.9 (dotted line). Using Eq.(^, it is easy to show that if we take 
Do = 3, for Dpi = 3 and Dpi -^ +oo, we have log j^ > — cx) and log p- = 0, respectively. 

FIG. 5. The generic shape of the potential Vii{D) for k = +1. The kinetic energy, D^, is positive 
in the classical region, Djj„i < D < Dt- The other values of D are corresponding to non-classical 
region in which the kinetic energy is negative. 

FIG. 6. The potential U{a) is shown by a solid line. The wave function of Hartle-Hawking 
(dashed curve) and that of Vilenkin (dotted curve) are shown. The real and imaginary parts of 
Vilenkin tunneling wave function are so indicated. The Hartle-Hawking wave function is real. 

FIG. 7. The generic shape of the potential f//(a). The growing wave function ^ i- (dashed line) 
and the decreasing wave function ^ i+ (dotted line) are shown. The under-barrier wave function, in 



the region 5 < a < ar is the sum of ^/- and ^1*/+ (see Eq jl25| ) . In the region ar < a < aim there is 
an oscillating wave function. For a > aim the decreasing wave function ^/_ is also indicated. This 
shape is drawn approximately and it shows only the behavior of the potential Ui{a). The similar 
shape can be drawn for the potential Uii{a) and its corresponding wave function, except for aum 
instead of aim- 
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